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ABSTRACT 

At a projected distance of ^26 pc from Sgr A*, the Arches cluster provides insight to star formation 
in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such 
as the cluster’s internal structure and orbital history are not well known. We present an astrometric 
and photometric study of the outer region of the Arches cluster (R > 6.25”) using HST WEC3IR. 
Using proper motions we calculate membership probabilities for stars down to E153M = 20 mag (^2.5 
Mq) over a 120” x 120” field of view, an area 144 times larger than previous astrometric studies of 
the cluster. We construct the radial profile of the Arches to a radius of 75” (^3 pc at 8 kpc), which 
can be well described by a single power law. Erom this profile we place a 3cr lower limit of 2.8 pc on 
the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of 
mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the 
orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced 
its closest approach to the GC between ^0.2 - 1 Myr ago. If accurate, this constraint indicates that 
the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, 
further simulations of clusters in the GC potential are required to interpret the observed profile with 
more confidence. 


1. INTRODUCTION 


The Arches cluster is a young (2-4 Myr; Najarro et al. 
2004, Martins et al. 2008) massive (^4-6 x 10^ Mq; 
Clarkson et al. 2012) star cluster near the center of the 
Milky Way. It has a projected distance of just ^26 pc 
from the supermassive black hole (SMBH) and is one 
of the most centrally concentrated star clusters in the 
Galaxy. Old enough to be free of its natal gas cloud 
and yet young enough to sample the full stellar mass 
range, the Arches cluster provides a unique opportu¬ 
nity to probe star formation and cluster evolution in the 
extreme Galactic Center (GC) environment. However, 
studies of the cluster are complicated by stellar crowd¬ 
ing and the high level of extinction which varies signifi¬ 
cantly acro ss the field (1.6 mag < A^s < 3.3 mag; Habibi 
et al.Jl2013). This effectively smears out the photometric 


properties of the cluster population, making it difficult to 
separate cluster members from field stars through pho¬ 
tometry alone. As a result, important questions about 
the cluster’s structure, initial mass function, and orbital 
history remain. 

The Arches is one of the closest examples of a young 
massive cluster (YMC) in a strong tidal field. Such ob¬ 
jects are not predicted to have long lifetimes, as simu¬ 
lated clusters near the GC sho w complete tidal disrup- 
tion on the order of ^10 Myr (Kim et al. 1999, 2000) 
or shorter when interactions with giant molecular clouds 
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are considered (|Kruijssen et al.|2014|) . Since the effects of 
tidal perturbations are mo st significant fm stars on the 
outskirts o f their clusters (Gnedin et al. 1999 Kiipper 
eFalllMQl ), measuring the structure of the outer region 
of the Arcnes offers insight to its past interactions. Eor 
example, it was long thought that the observed tidal ra¬ 
dius (i.e., limiting radius) of a cluster, where the stellar 
density drops to zero, should correspond to its theoret¬ 
ical tidal radius (i.e., Jacobi radius), where the gravita¬ 
tional acceleration of the clu ster equals the tida l accel¬ 
eration of its parent galaxy (|von Hoerner 1957). Early 
Eokker-Planck simulations of clusters on eccentric orbits 
further suggested that the tidal radius imposed by the 
strong est tides at perig alacticon should persist to later 
times (Oh & Lin 1992). However, more recent N-body 
simulations show that perigalacticon passage does not 
cleanly truncate a cluster, but rather results in an ex¬ 
tended radial profile that approaches a power law (Oh 


et al. 11199^ |Jc)hnston et al.||1999| |Peharrubia et al.||26u9[ 


.poiOD 

le slope 


Kiipper et al.||2010|). The “extratidal” stars can have a 


ditterent profile slope than the rest of the cluster, creating 
a break in the profile at a radius which may be related to 


the t i me since perigalac ticon passage (jPeharrubia et al. 
2009 Lokas et al.||2Q13 ) 


Measurements of globular cluster profiles have revealed 
the presence of ext ratidal stars, often associated with 
tidal tail structures dOdenkirchen et al.|2001||Siegel et al. 


20011 lOdenkirchen et aL|20031|HeloknrQv et aL|2006[|bol- 


li ma et al.|2011||Chun e t al.| 2015). However, as di s cussed 
bylKiipper et al. (|2010|) ana Garballo-Hello et al. (2012), 
even clusters that are not exposed to varying external 
tidal fields can also form extratidal structures through 
two-body relaxation, as stars which become energetically 
unbound from the cluster do not escape instantaneously 
but rather on a time-scale which is dependent on their or 


bital parameters 

(lEukushige & Heggie||2000 Baumgardt 

& Makino||2003 

Zotos||2015p. Both of these mechanisms 
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affect the radial profile of globular clusters, which have 
ages larger than their relaxation times. On the other 
hand, two-body relaxation is negligible for YMCs that 
are much younger than than their relaxation times, and 
so the impact of tidal perturbations on their profiles 
should be easier to isolate. Given the high likelihood of 
strong tidal interactions with the GC, the Arches offers 
a promising opportunity to measure such a feature. In 
addition, a detailed understanding of the Arches cluster 
profile is necessary for assessing the impact of dynami¬ 
cal effects on the present-day mass function (i.e. mass 
segregation and tidal stripping). 

Insight into the past tidal interactions of the Arches 
cluster also provides valuable inform ation about its or¬ 
bit. While its bulk proper motion ( Stolte et al. 200^ 


Clark son et al. 2012) and doppler velocity (f'iger et al 
2U02[ ) have been measured, the line-of-sight distance is 


unknown, preventing a unique orbital solution. As a re¬ 
sult, the birth environment of the cluster and its rela¬ 
tion to the nearby Quintuplet cluster are not well under¬ 
stood. For example, it has been suggested that both the 
Arches and Quintuplet formed in a collision between gas 
clouds along the Xi and X 2 orbit families in the Galac- 



to the GC. It has also been suggested that the Arches 
may be a possible source of the isolated massive stars 
observed near the GC, depending on its orbit and how 
much tidal strippin g has occurred (jMauerhan et al.|2010 
Habibi et al.||2Q14 ). If the time since the last tidal per 


turbation (presumably occurring at the closest approach 
to the GC) can be establishe d, the set of possibl e orbits 
for the Arches calculated by Stolte et al. (2008) can be 
significantly restricted. 

We present the radial profile of the Arches cluster out 
to large cluster radii (^75”, or ^3 pc at 8 kp(Q and in¬ 
vestigate the structure of the outer region of tne cluster 
for the first time. We use stellar proper motions rather 
than photometry to calculate cluster membership prob¬ 
abilities, avoiding many of the difficulties introduced by 
differential reddening. The effectiveness of this method 
on the Arch es cluster was demonst rated by [Stolte et al.| 
(2008) and Clarkson et al. (2012), who used ground- 
based adaptive optics (AO) observations to measure the 
cluster’s bulk proper motion and identify members in the 
central 10” x 10” region. In this paper, we conduct an 
astrometric study of the Arches cluster using the Hubble 
Space Telescope (HST), which provides high astrometric 
precision over a field of view 144 times larger than these 
previous studies. This allows us to measure the radial 
profile to beyond the predicted tidal radiu s (25” - 60”; 
Kim et al. 2000 Portegies Zwart et al. 2002|). In addition. 


we examine the degree of mass segregation throughout 
the cluster and search for the presence of tidal tails. The 
consequences of our results are discussed in relation to 
the orbital history of the Arches cluster. 


^ All distances throughout this paper assume a distance of 8 kpc 
to the Arches cluster. 



Figure 1. Three color image of the Arches Cluster, with 
F127M = blue, F139M = green, and F153M = red. Significant 
differential extinction is apparent from the changing density of the 
field stars. The hole in the lower left side of the image is due to 
a known defect in the WFC3IR chip. The proper motion of the 
cluster (labeled with green arrow) is very nearly parallel to the 
Galactic plane. 


2. METHODS 

2.1. Observations and Measurements 

We observed the Arches cluster with the Hubble Space 
Telescope (HST) WFC3IR camera using the F127M, 
F139M, and F153M filters (1.27 /im, 1.39 /im, and 1.53 
/im, respectively; PI: Ghez, ID: 11671, 12318, 12667^ 
A summary of the observations is provided in Table [l] 
These observations have a field of view of 120” x 120” 
and are centered at (a(J2000) = 266.4604, (5^2000) = - 
28.8222 with a position angle of -45° (Figure [^. Astrom¬ 
etry is performed on the F153M observations, which were 
obtained in three epochs over a two year baseline between 
2010 - 2012. Of this filter set, F153M was chosen for as¬ 
trometry because it provides the optimal combination of 
limited saturation and a well-sampled point spread func¬ 
tion (PSF) with a FWHM ^0.17” (1.4 pix, scale = 0.121” 
pix“^). High astrometric and photometric precision is 
achieved by observing at the same position angle and 
pixel position across epochs and using a dense, sub-pixel 
spiral dithering pattern within each epoch. F127M and 
F139M observations were only obtained in 2010 using a 
simpler dither pattern for the purpose of color informa¬ 
tion to derive extinction. 

Each frame is processed using the standard HST 
pipeline, which produces a fit image which has been 
flat-fielded and bias-subtracted. We use a combination 
of public and custom software to extract high-precision 
astrometry and photometry via PSF-fitting with a vari¬ 
able PSF model. While this method has been previously 
implemented on HST optical observations of Globular 
Clusters, this is the first time it has been applied to 
WFC3IR observations. We produce a list of individual 
stellar measurements for each filter/epoch using the pro¬ 
gram KS2^ a generalization of the software devel oped to 
reduce the Globular-Cluster Treasury Program (Ander- 
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Table 1 

HST WFC3IR Observations 


Date 

Filter 

^images 

Total Exp. Time (s) 

Depth (mag)^ 

Pos Error (mas)® 


Phot Error (mag)^ 

2010.6150 

F127M 

12 

7200 

23.63 

0.60 


0.02 1 

2010.6148 

F13yM 

10 

3500 

23.29 

0.90 


0.03 - 

2010.6043 

Flb3M 

21 

V350 

23.31 

0.88 


0.05 1 


2012.6156" 

F153M 

21 

V350 

23.31 - 

0.88 


0.05 1 

^ Median value at mag 

= 20 in n 

espective filter 





Estimated from 95th percentile of extracted stellar magnitudes. | 


son et al. 2008). Star positions are transformed to an 
arbitrary astrometric reference frame where the net mo¬ 
tion of the cluster plus field is 0 mas yr“^ using general 
6-parameter linear transformations that can be described 
as a 2D translation, rotation, plate scale, and shear for 
each image. Approximately ^50,000 stars are measured 
in each filter/epoch. The photometry is calibrated to the 
standard Vega magnitude system using the zero points 
derived for the WFC3IR camersj^ A detailed description 
of the data reduction and measurement process, as well 
as an analysis of the astrometric and photometric errors, 
is provided in Appendix [A| 

We calculate proper motions for the stars that are 
detected in all three F153M epochs, using linear fits 
weighted by the individual astrometric errors. Proper 
motion uncertainty as a function of observed F153M 
magnitude is presented in Figure Several tests were 
conducted to confirm the validity of these errors (see 
Appendix 0 . Previous studies of the bulk proper mo¬ 
tion of the Arches c l uster relati ve to the field popnla- 
tion by [Stolte et al. (2008) and Clarkson et al. (2012) 
revealed that a precision of ^0.8 mas yr“^ is needed 
in order to reliably separate cluster members from field 
stars. As a conservative error cut, we restrict the forth¬ 
coming analysis to stars with a proper motion precision 
of 0.65 mas yr“^ or better, which we achieve down to 
F153M « 20 mag. At the average distance and redden¬ 
ing of the Arches cluster this corresponds to roughly 2.5 
Mq. We additionally require a minimum photometric 
precision of 0.06 mags in each F153M epoch to ensure 
high-quality results. We measure ^26,000 proper mo¬ 
tions, ^6000 of which pass these error cuts. The kine¬ 
matic distinction between cluster and field stars is clearly 
seen in a vector point diagram (VPD; Figure [^. Note 
that the proper motions have been rotated from image 
coordinates into projected equatorial coordinates where 
the two-dimensional proper motion vector of any star is 
fi = [/ic^cos(5, jj^s]. Proper motions are also shifted into 
a reference frame where the cluster is at rest, estimated 
from the mean motion of stars within the central 10”xl0” 
region of the cluster. This sample is an order of magni¬ 


tude l arger than the sample analyzed by [Clarkson et al. 
(2012), who had a much smaller field of view. 


2.2. Cluster and Field Populations 

In order to calculate cluster membership probabilities 
we must first characterize the kinematic distributions of 
the cluster and field star populations. Previous studies of 
the Arches cluster have assumed that the field kinematics 
can be modeled as a single elliptical Gaussian distribu- 


® As of August 2014; http://www.stsci.edu/hst/wfc3/phot_zp_lbn 


10 ^ 



Observed F153M Mag 


Figure 2. Proper Motion error vs. F153M magnitude. The red 
line marks our proper motion error cut of 0.65 mas yr“^. Only 
stars below this cut are included in our analysis. The green line 
denotes F153M = 20 mag, which corresponds to ~2.5 Mq at the 
approximate distance and average reddening of the Arches cluster. 


tion (Clarkson et al.|2012). The field kinematic distribu¬ 
tion is elliptical because it is primarily composed of stars 
in the Galactic bulge that exhibit a larger velocity disper¬ 
sion along the Galactic plane than perpendicular to it, a 


consequence of coherent rotation (Clarkson et al. 2008 



jy a single 

account for this complexity, we ad opt a Normal Mixture 
Model (McLaughlin & Peel 2000) to simultaneously fit 
multiple Gaussians to the observed VPD. A more com¬ 
plete description of the field requires modeling the stellar 
density, kinematics, and reddening at all distances along 
our line of sight towards the Arches. Such analysis would 
be valuable for exploring Galactic structure but is beyond 
the scope of this paper. 

We construct a likelihood function for each star in the 
sample from from the sum of K Gaussian components: 
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Figure 3. Vector point diagram of the ~6000 stars included in our 
analysis. Proper motions are in the reference frame of the cluster. 
Cluster members appear as a distinct clump of stars at {/j^aCOsS, 
11 s) = (0, 0), while the field stars are spread along the Galactic 
plane. 


^ ^ 1 / 1 

( 1 ) 

N 

C = X{L{^li) 


where fii is the proper motion of the ith star, tt/c is 
the fraction of total stars in the kth Gaussian such that 
X]/c=o = I 5 M/c is the velocity centroid of the kth Gaus¬ 
sian, and Tij^i is the covariance of the kth Gaussian and 
ith star. The t otal likelihood over th e sample of N stars 
is jC. Following Clarkson et ah ( 2012 ), we add the covari¬ 
ance matrices of the population model and stellar proper 
motion uncertainties such that where Si 

is the velocity error matrix (assumed to be diagonal with 
velocity error components and cr^^) and Zj^ is the 

covariance matrix of the kth Gaussian fit. 

With the likelihood function defined, we can determine 
the global kinematic parameters of the cluster and field 
populations through Bayesian inference using Bayes’ the¬ 


orem: 


P{tt,h,Z\h,S) = 


Pin, g|7r, /Lt, Z)P{tt, h, Z) 
P{ti,S) 


( 2 ) 


where P( 7 r, 7 x, Z\ii^ S) is the posterior probability of our 
model parameters tt, the set of nk values; the set of 
Gaussian velocity centroids; and Z, the set of Gaussian 
covariance matrices given the observed stellar velocities 
fi and velocity error matrix S. S\7V,JI, Z) is the 


probability of the observed stellar velocity distribution 
given the model, and P( 7 r, 7 I, Z) is the prior probability 
of the model. In this case, S\7V^JI, Z) is the total 
likelihood jC defined in Equation 

To find the posterior probability distribution we use 
Multinesk a publicly available nested sampling algorithm 
which serves as an alternative to Markov Ghain Monte 
Carlo (MCMC) al gorithms when ex ploring multi-modal 
parameter spaces ( [Feroz et al.|2QQ9 ). This iterative tech¬ 
nique calculates the posterior probability at a fixed num¬ 
ber of points in the parameter space and identifies possi¬ 
ble peaks, restricting subsequent sampling to the regions 
around these peaks until the change in evidence drops 
below a user-defined tolerance level. Multiple peaks can 
be identified and evaluated, resulting in increased sam¬ 
pling efficiency with complicated parameter spaces. We 
run the algorithm us ing the python module PyMultinest 
(jBuchner et al.||2014|). 

We hud that the cluster and field populations can be 
well described with a 4-Gaussian mixture model, with 
one Gaussian describing the cluster and the other three 
describing the field (Figure 0. The use of this model is 
justifi ed by the Bayesian Information Criterion ( [Schwarz' 


1978), which is minimized compared to less coniplicated 


(3-Gaussian) or more complicated (5-Gaussian) mixture 
models. We require the clus ter Gaussian to be cir cular, 
consistent with the results of Clarkson et al. ( 2012 ), and 
adopt a prior to roughly constrain its location around 
(va;, Yy) = (0,0). The parameters for the field Gaussians 
as well as the remaining parameters for the cluster Gaus¬ 
sian are unconstrained. The one-dimensional posterior 
distributions are well described by a Gaussian function, 
Iwhich is used to determine the best-fit value and error 
for that parameter. A summary of the parameter priors, 

best-fit results, and errors is provided in Table _ 

These res ults can be compared with those of fClarksonj 
et al. ( 2012 ), who examine the kinematics of the inner 10 ” 
X 10 ” of the cluster using ground-based AO observations. 
Their measurements have a higher precision but much 
smaller field of view than our observations. We obtain 
a velocity dispersion of 0.18 ± 0.02 mas yr“^ for the 
cluster, which is con sistent with the measu rement of 0.15 
± 0.01 mas yr“^ by Clarkson et al. ( 2012 ) within errors. 
This agreement comes despite using fully independent 
data sets which focus on different regions of the cluster. 

However, there is less agreement on the bul k motion 
of the Arch es relative to the field population. [Clarkson 


et al. (2012) model the field using a single elliptical Gaus- 


sian distribution with a velocity center offset by 4.39 ± 
0.38 mas yr“^ from the cluster. In this study we model 
the field using 3 elliptical Gaussians, and our fits indicate 
tha t all have smaller mot ions relative to the cluster than 
the[Clarkson et al.[([ 2012 [) result (see Table [^. If we cal- 
culate the average motion of the field from the average 
sum of the 3 field Gaussians, then we get an overall field 
motion of 2.83 ± 0.33 mas yr“^ relative to the cluster, 
which is substantially and significantly lower than the 


Clarkson et al. ([ 2012 ) result. We discuss this discrep- 


ancy further in §4.1 [ 

2.3. Extinction Map Using Red Clump Stars 

Taking advantage of the high photometric precision 
of HST, we use red clump (RC) stars in the Galactic 
bulge to measure the extinction across the Arches clus- 
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Table 2 

Cluster and Field Population Model: Free Parameters, Priors, and Results 


I Cluster Caussian Field Gaussian 1 Field Gaussian 2 Field Caussian 3 

Parameter^ Prior^ Result Prior Result Prior Result Prior Result 



'^k 

U(0, 1) 

0.08 


0.01 

U(0, 1) 

0.25 


0.03 

U(0, 1) 

0.42 


0.04 

U(0, 

1) 

0.25 


0.03 

Ta.,k 

(mas yr“^) 

G(0, 0.2) 

0.06 


0.03 

U(-4, 12) 

-0.26 


0.06 

U(-4, 12) 

-1.19 


0.11 

U(-4, 

12) 

-1.15 


0.15 

TS,k 

(mas yr“^) 

G(0, 0.2) 

0.06 


0.02 

U(-4, 12) 

-0.77 


0.09 

U(-4, 12) 

-3.40 


0.25 

U(-4, 

12) 

-3.05 


0.20 

^a,k 

(mas yr“^) 

U(0, 8) 

0.18 


0.02 

U(0, 8) 

1.27 


0.09 

U(0, 8) 

2.73 


0.12 

U(0, 

8) 

3.20 


0.13 

^b,k 

(mas yr“^) 


0.18 


0.02 

U(0, 4) 

0.60 


0.05 

U(0, 4) 

1.30 


0.08 

U(0, 

4) 

3.05 


0.13 

Ok (rad) 

— 


0 


U(0, tt) 

1.16 


0.05 

U(0, tt) 

1.20 


0.03 

U(0, 

tt) 


0 



^ Description of parameters: = fraction of stars in Gaussian; = RA-velocity centroid of Gaussian; = DEC- 

velocity centroid of Gaussian; aa^k = semi-major axis of Gaussian; aij^k = semi-minor axis of Gaussian; 9k = angle between 
cFa^k s-iid the RA-axis 

^ Uniform distributions: U(min, max), where min and max are bounds of the distribution; Gaussian distributions: G(yLi, cr), 
where fi is the mean and a is the standard deviation 





Figure 4. Top: Vector point diagram of our sample with the fitted 1- and 2cr distributions of the cluster Gaussian (red) and field 
Gaussians (blue, green, cyan corresponding to field Gaussians 1, 2, and 3 in Tablerespectively). The left plot shows all stars in the 
field, demonstrating the extension of the field populations in the direction of the Galactic plane. The right plot is a zoomed-in view of the 
cluster population, readily apparent as a tight clump of stars moving with a common motion relative to the field. Bottom: Proper motion 
distribution of the stars in projected equatorial coordinates {left: RA, right: Dec). The predicted distribution of the Normal Mixture 
Model (red) is found to be a good match to the observed stars (black). 
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ter field. This provides an alternative to “sliding” ap¬ 
parent cluster members along their reddening vector in 
a color-magnitude diagram (CMD) to a theoretical clus¬ 
ter isochrone to measure extinction, as has been done in 
previous studies of the Arches cluster (Kim et ah 2006 


Espinoza et al.j2009t|Habibi et al.|j2013|). This (JMJJ slid¬ 

ing method is prone to held contamination and isochrone 
uncertainties, especially for the pre-main sequence at low 
to intermediate masses. On the other hand, stellar evo¬ 
lution theory and observations show that RC stars ex¬ 
hibit well-defined luminosities and colors which do not 


vary significantly with age or metallicity (Castellani et al. 


1992 Paczynski & Stanek 1998 Stanek et al.|2000p, mak¬ 
ing them usetul calibrators to measure extinction. This 
is especially true near the GC, where the relatively high 
density of RC stars in the bulge population makes i t pos¬ 
sible to create reddening m aps of different regions (Sumi 
2004| |Schodel et al.||2010|) . Though the line-of-sight po¬ 
sition oTtEe'TfrcEes'^witirrespect to the bulge RC stars is 
uncertain, we assume that all of the extinction is caused 
by foreground material and so the RC population ex¬ 
hibits similar extinction as the cluster members them¬ 
selves. With this approach we create the first RC-based 
extinction map of the Arches cluster. 

Bulge RC stars are readily identified as a narrow pop¬ 
ulation spread along the reddening vector in the F127M 
vs. F127M - F153M color-magnitude diagram (Figure 
1^. The spread of this population is primarily caused 
oy differential reddening, which smears out what would 
normally be a tight clump of stars. In order to isolate 
RC stars, w e use a PHOENIX model atmosphere (Allard 
et al.|2011 ) with typical RC para meters at solar metallic- 
/I'e// = 4700 K, log g = 2.40; |Mishenina et al.|2006) 
to calculate the F127 vs. F1271vl - 4’I53M reddening 
vector using the GC extinction law of |Nishiyama et ah 
(2009). Keeping the slope of the reddening vector hxed, 
the y-intercept is fit to stars which fall within a broad 
area in the CMD around the RC population. We identify 
RC stars as those within a rectangle with the long axis 
centered on the reddening vector with the length of the 
short axis defined by the least crowded section of the RC 
bar (F127M - F153M « 2.7). This corresponds to a width 
of constant val ue A F127M = 0.7 mag. Identified clus¬ 
ter members (§ 3.1) are removed from this sample. The 
extinction of each identified RC star is taken from the 
nearest point on the reddening vector in color-magnitude 
space. 

We measure the extinction for 1027 RC stars identified 
across the field. These values are spatially interpolated 
using a 5th order bivariate spline to map the extinction 
at every position (Figure [^. T he typical error is (TAks 
= 0.10 mags, as derived i]:i §3.1| Extinction values range 
from 1.8 < Aks < 3.0 with a median of Aks = 2.4 for 
cluster members. This range is in aCT eement with the 
reddening map of Ha bibi et al. (|2013|) , who find 1.6 < 
Aks < 3.3 also using a|JNishiyama et al. (2009) reddening 
law. 


2.4. Completeness Analysis 

In order to accurately measure the radial density pro¬ 
file of the Arches Cluster we must conduct an extensive 
completeness analysis on our astrometry pipeline. In ad¬ 
dition to the sensitivity threshold of our observations, 
stars may be missed due to source confusion and prox¬ 


imity to bright and/or saturated stars. These effects are 
especially relevant for the dense central region of the clus¬ 
ter. To quantify our completeness we perform an arti¬ 
ficial star injection and recovery test, planting 400,000 
stars in each image and determining which are recov¬ 
ered to sufficient accuracy and precision as a function of 
spatial position and magnitude. The magnitudes of the 
artificial stars are drawn from the observed CMD of the 
field in order to best simulate the photometric properties 
of the observed stars. These magnitudes are then per¬ 
turbed by a random amount reflecting the photometric 
uncertainty (assumed to be Gaussian distributed) of the 
real star they are simulating. The same set of artificial 
stars is applied to all observations. This analysis assumes 
that the artificial star measurement errors match those 
of the observed stars, which we test in Appendix [B] 

The conditions an artificial star must fulfill in order 
to be considered as recovered matches the criteria ap¬ 
plied to the real data. Within a given epoch, a recovered 
artificial star must: 1) be detected in at least 75% of 
the images within that epoch; 2) have position and mag¬ 
nitude errors less than 1.5 mas (required for a proper 
motion precision better than or equal to 0.65 mas yr“^) 
and 0.06 mag, respectively; and 3) have a measured po¬ 
sition and magnitude within 0.5 pix (60 mas) and 0.5 
mag of the planted values to guard against misidentifi- 
cation. In addition, artificial stars must be recovered in 
all 3 F153M epochs, which is required of the obs erved 
stars in order to derive their proper motions (§ |2.1[ ). Af¬ 
ter detection/non-detection, the extinction map is used 
to differentially de-redden the artificial stars to the mean 
extinction of the cluster {Aks = 2.4 mag). The fraction 
of recovered artificial stars to the total number of planted 
artificial stars represents the completeness fraction as a 
function of position and differentially de-reddened mag¬ 
nitude. 

The resulting completeness curves as a function of dif¬ 
ferentially de-reddened magnitude and of radius are pre¬ 
sented in Figure Over the full field we achieve greater 
than 50% completeness down to F153M = 20 mag. How¬ 
ever, the completeness in the inner 6.25” (^0.25 pc) of 
the cluster is significantly lower due to stellar crowding, 
falling to 30% by F153M = 18.5 mag. As a result, we 
restrict the following analysis to observed stars with R > 
6.25” and differentially de-reddened magnitudes brighter 
than F153M = 20 mag. 


3. RESULTS 

3.1. Cluster Membership 

With the kinematic pro perti es of the cluster and field 
populations determined (§ 2.2), we calculate the proba¬ 
bility of cluster membership tor each star based on its 
proper motion: 


pi 

member 


oPc + Ef T^kP], 


(3) 


where tTc and tt/c are the fraction of total stars in the 
cluster and kth field Gaussian, respectively, and PJ, and 
are the probability of ith star being part of the clus¬ 
ter and kth field Gaussian, respectively. A histogram of 
the resulting cluster membership probabilities is shown 
in Figure In the following analysis we include all 
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Figure 5. The identification of Red Clump (RC) stars in the F127M vs. F127M - F153M color-magnitude diagram, as described in § |2.3| 
Left: The full field CMD, with the stars used in the initial fit of the RC reddening vector in red and the reddening vector itself in blue. 
Right: A zoomed-in view of the RC population, with the RC reddening vector and final identification criterion shown with blue and green 
lines, respectively. Stars falling between the green lines (red points) are identified as RC stars and are used to make the extinction map. 



X position Offset (arcsec) 


Figure 6. The extinction map created via the spatial interpola¬ 
tion of the Red Clump extinction values as de scribed in the text. 
The IR reddening law of|Nishiyama et al.|(|2009|) is used to calculate 
Ax at different wavelengths. Fositions are given with respect to the 
cluster center located at (x, y) = (0, 0), and the axes are oriented 
same manner as Figure [R Strong cluster candidates (Pmember > 
0.7) have an average redaening of A^s = 2.4 mag. 


stars with Vmember > 0-3, Weighted by their individual 
membership probabilities. This criteria selects 701 stars 
which represent 446.8 “cluster members” based on the 
sum of the cluster membership probabilities. We con¬ 
sider stars with Pmember > 0-7 as strong cluster candi¬ 
dates, whose distributions in position and velocity space 


are shown in Figure 

By applying the extinction map derived in §2.3| to the 
strong cluster candidates we create a differentially de- 
reddened F127M vs F127M - F153M CMD of the clus¬ 
ter (Figure [Tq|) . The improvement relative to the un¬ 
corrected CMD is noticeable in both the overall color 
dispersion and definition of the blue edge. Stars with 
colors more blue than the blue edge are very likely field 
contaminants, while the the scatter along the redward 
edge of the cluster may be caused by intrinsic reddening 
of the objects themselve s, perhaps due to circumstellar 
disks (Stolte et al.ll 


The remaining color dispersion of the differentially de- 
reddened CMD provides an estimate of the uncertainties 
of the extinction map. Between 16 < F127M <21 mag 
the median color dispersion is 0.36 mags. Assuming that 
the photometric uncertainties in each filt er are negligible 
and adopting the Nishiyama et al. (2009|) extinction law, 
this color dispersion corresponds to an extinction error 
of (7Afi53m = 0-13 mags {(Taks = 0.10 mags). 

The mean extinction of the strong cluster candidates is 
Krs = 2.42 ± 0.14 mag. The inner region of the cluster 
(R < 0.4 pc) exhibits a tight range of reddening values 
from 2.33 < Ars < 2.53 mag, while the outer region 
of the cluster (R > 0.4 pc) exhibits a much wider range 
from 2.04 < Krs <2.76 mag. In the literature, there are 
variations in the measured extinction for the Arches due 
to differ ent methodologies and reddening laws. Using 
the same Nishiyama et al. (|2009) ext i nction law with the 
CMD sliding method, Habibi et al.|(2013) find A^s = 
2.6 ± 0.2 mag for 0.2 pc < R < 0.4 pc and Ars = 2.6 
± 0.3 mag for 0. 4 pc< R < 1.5 pc. Also usin g the CMD 


sliding method, Espinoza et al. (2009) and Kim et al 
(|2006]) find higher values of A/y.s = 2.97 and 3.1 for R 
< 0.4 pc using the extinction laws of Fitzpatrick (2004) 
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Figure 7. Completeness as a function of observed magnitude and radius. Left: Completeness as a function of differentially de-reddened 
magnitude {A^s = 2.4 mag) for the F127M, F139M, and F153M filters. These completeness values are calculated over the entire field. 
Right: The F153M completeness as a function of differentially de-reddened magnitude in different radius bins. The blue line marks a 
minimum completeness of 30%, which is achieved down to F153M = 20 mag for R > 0.25 pc. This sets the faint-end magnitude limit and 
inner radius limit for the radial profile. 





1000 


8.0 0.2 0.4 0.6 0.8 1.0 

Cluster Membership Probability 


Figure 8. A histogram of the membership probabilities obtained 
for the sample. All objects with P cluster ^ 0-3 (red line, inset plot) 
are considered in the profile analysis, weighted by their membership 
probability. Of ~6000 stars examined, 701 meet this criterion, with 
membership probabilities that sum to 446.8. 


and Rieke et al. (1989), respectively. Our measurements 
are consistent with previous measurements made using 
the same extinction law. 


3.2. Radial Density Profile 

Using stars with Prnember > 0-3, we construct the ra¬ 
dial profile of the Arches Cluster us i ng th e Bayesian 
methodology described by |Do et al.| (|2013|). This al¬ 
lows us to construct an un-binned profile that simultane¬ 
ously incorporates cluster membership probabilities, im¬ 
age completeness, and geometric area corrections at large 
radii to a ccou nt for incomplete area coverage. As dis¬ 
cussed in §2.4[ only stars with differentially de-reddened 


magnitudes of F153M < 20 mag and R > 6.25” (0.25 
pc) are considered. We adopt a single power law as our 
likelihood function: 

Li{r,r,b)=Aorr^ + b (4) 

where is the radius of the ith star and the field con¬ 
tamination b is assumed to be constant across the im¬ 
age. The profile amplitude Aq is calculated such that 
the integral of the radial profile yields the total number 
of cluster members observed after membership probabil¬ 
ity, completeness, and area corrections. 

The total likelihood jC is then the product of the indi¬ 
vidual likelihoods for N total stars: 

N 

£ = J2Mr)iogLi{r,r,b) (5) 


Wi{r) 


Pi 

Aiir)Ciir) 


where Pi is the membership probability of the ith star, 
Afii) is the relative fraction of observed area at the star’s 
radius relative to an infinite field of view {Ai = 1.0 for 
0 < r < 60”, Ai < 1.0 for r > 60”), and Cfii) is the 
completeness at that star’s radius. A summary of our 
best-fit model and subsequent results are presented in 
Table and Figure [m A binned profile is included for 
comparison, with errors calculated from the Poisson un¬ 
certainties in the completeness correction and observed 
profile as well as the uncertainties in the extinction map. 
These are captured by recalculating the stellar density in 
each radius bin using a magnitude cut brighter or fainter 
than F153M = 20 mag by the map error value {cfafi^sm 
= 0.18 mag). The half-light radius of the profile i s (3.48 
pc, lar gely consistent with previous studies (0.4 pc,|Figer 


et al.| 1999). The bivariate posterior distributions for 
these parameters are presented in Appendix \C\ 

Included in the right panel of Figure is the radial 
profile for t he inner part of the Arches n*om [Espinoz^ 


et al. (2009), which spans a stellar mass range of 10 Mq 
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Figure 9. The spatial (left) and kinematic (right) positions of strong cluster candidates (Prnember > 0-7), marked as red points, compared 
to the rest of the sample in black points. The spatial positions are plotted in arcseconds relative to the cluster center and are in image 
coordinates (same orientation as Figure^. Proper motions are plotted in projected equatorial coordinates. The radial profile of the Arches 
includes all stars with Prnember > 0-3, larger sample than is shown here. 



F127M-F153M F127M - F153M F127M-F153M 


Figure 10. The F127M vs. F127M - F153M color-magnitude diagram for the full sample (/e/t), strong cluster candidates (Pmember > 0-7, 
middle), and strong cluster candidates after being differentially de-reddened using the extinction map (right). The reddening correction 
noticeably tightens the color dispersion and blue edge of the population. The differentially de-reddened CMD is consistent with a theoretical 
2.5 Myr cl uster isochrone at 8000 pc with Aks = 2.4 mag, overlaid in red. The isochrone is c reated using the pre-main s equence evolutionary 
models of|Siess et al.| |2000| for M < 7 Mq and the main sequence models with rotation of |Meynet & Maeder| (|2003|) for M > 9 Mq , with 
an interpolation between the models over the missing mass range. 


< M < 120 Mq out to R = 0.4 pc. There is good agree¬ 
ment between the sh ape of the two profiles , though the 
absolute values of the Espinoza et al. (2009) profile must 
be scaled. This is necessary due to differences in sensitiv¬ 
ity and treatment of cluster membership. We note that 
our profile spans far beyond the limits of previous astro- 
metric studies of the Arches cluster, which are restricted 
to R < 0.2 pc. We leave the combination of these astro- 
metric data sets and the presented data set to a future 
paper. 

To quantitatively assess whether the power-law model 


is an appropriate one for the observed profile, we con¬ 
duct a p osterior predictive an alysis using as the test 
statistic ( Gelman et al.|[2M3 ). We randomly select 1000 
sets of model parameters from the joint posterior distri¬ 
bution and generate artificial binned profiles from these 
models. Each data point within the artificial profiles is 
shifted by an offset randomly drawn from a normal dis¬ 
tribution with a width equal to the uncertainty in that 
value, determined from the combination of the poisson 
uncertainty and the uncertainty in the completeness cor¬ 
rection. We then calculate a value for each binned 
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Table 3 

Power-Law Profile Fit Results 



Bin 


Power-law slope 

F 

Field Contamination 
b (stars / pc^) 

Normalization Constant 
Aq (stars / pc^) 

Full Cluster 


451.0 

2.06 ± 0.17 

2.52 ± 1.32 

23.09 ± 3.5 

Split by Mass (2-bin) 

High Mass^ 

106.5 

2.70 ± 0.35 

0.64 ± 0.48 

3.50 ± 1.22 


Low Mass‘d 

354.5 

1.75 ± 0.15 

0.78 ± 1.37 

20.33 ± 2.75 

Split by Mass (3-bin) 

High Mass^ 

129.6 

2.75 ± 0.37 

1.47 ± 0.48 

3.36 ± 1.26 

Intermediate Mass® 

163.6 

2.00 ± 0.28 

2.04 ± 0.56 

6.38 ± 1.46 


Low Mass^ 

165.7 

2.29 ± 0.30 

2.09 ± 0.53 

6.47 ± 1.48 

Split by Direction^: 

Parallel 

225.7 

1.86 ± 0.17 

2.32 ± 0.88 

22.57 ± 1.9 

Prior^ 

Perpendicular 

226.0 

2.19 ± 0.18 
U(0.5, 4.5) 

2.29 ± 0.82 
U(0,8), G(2.52,1.32)' 

21.35 ± 2.43 


^ Weighted by membership probability and corrected for completeness. 

F153M < 17 mag (M > ~13 M©) 

17 < F153M < 20.0 mag (~2.5 M© < M < ~13 M©) 

F153M < 17.3 mag (M > ~12 M©) 

® 17.3 < F153M < 18.8 mag (~6 M© < M < ~12 M©) 
f 18.8 < F153M < 20.0 mag (~2.5 M© < M < ~6 M©) 
s Relative to the direction of the Arches cluster orbit 

^ Uniform distributions: U(min, max), where min and max are bounds of the distribution; Gaussian distributions: G(yLi,(7), 
where /j, is the mean and a is the standard deviation 

^ Adopted U(0,8) for the full cluster profile fit, G(2.52, 1.32) for the directional profile fit 


profile with respect to the best-fit model to the observa¬ 
tions: 


x\T,b) = 


^pbin _ p'model'^2 


i=i 




(6) 


where is the jth point in the binned profile with un¬ 
certainty aj and is the value predicted for the jth 

bin by the best-fit model. We find that only 1% of these 
values are lower than the value for the observed 
binned profile and conclude that a power-law model is a 
good fit to the data. We emphasize that this x^ statistic 
and binned profiles are used to test the validity of the 
model, not to fit the model itself. 

Previous studies of other YMCs have shown that these 
objects often appear to have exte nded radial profiles 
without signs of tidal truncation (|Elson et al.| |1987 


0 . 2 , which is consistent with the range and median val¬ 
ues of 2.01 - 3.79 and 2.5 9 determined from a sample o f 
LMC and SMC YMCs by|Mackey & Gilmorej (|2003a|b|). 
Given that we have no information about Vc from our 
profile and the marginal difference between this profile 
and the power-law fit, we proceed with the single power- 
law model. 

3 . 2 . 1 . Mass Segregation 

Evidence for mass segregation in the Arches cluster has 
been found in the f lattening of the mass function toward 
the cluster center (|Figer et al. 111999111 iS tolte et al .||2002 
Espinoza et al.||200^ Habibi et al.||2013p and a shallower 
radial profile for stars between ^10 - 30 M^;^ compared 
to stars between ^30 - 120 Mq (Espinoza et al. 2009). 


Mackey fc Gi 


____ ilmore|2003a|b||McLaughT^^^ ^ ___ ^ 

2005D. These studies ht the prohles using a model dehned stellar crowding ( [AscensQ et al. [2009|) 


. & van der Mare. 


However, in addition to being dependent on photometric 
cluster membership, these results rely on measurements 
in the dense innermost regions of the cluster (R < 15”) 
where completen ess is lowest for low-m ass stars due to 

We avoid this 


By Elson et al. (1987), hereafter referred to as an EEE87 
profile model. We ht our profile with this model, adding 
a constant term b for field contamination: 


S(r) = Eo 1 


-7/2 


( 7 ) 


where a is related to the core radius of the cluster: 


(2^/1 _ 1) 


1/2 


( 8 ) 


We use the same Bayesian framework as described 
above, only with this profile as the likelihood function. 
Since incompleteness prevents our profile from stretch¬ 
ing into the core region o f the Arches we adopt the core 
radius Vc determined by Espinoza et al. (2009) of 0.14 
± 0.05 pc in our model. 4'he consequent ht is nearly 
identical to the single power law model fit, with param¬ 
eter values provided in Table We obtain 7 = 2.3 ± 


inner region and instead examine mass segregation in 
the less-den se outer regions of th e cluster. 

Following |Espinoza et al.|(|2009|), we separate our radial 
profile as a function of ditterentially de-reddened magni¬ 
tude to test for mass segregation. These magnitudes are 
a good proxy for stellar mass, as only a small fraction 
of stars in the Arches (^ 6 %) have been found to exhibit 
IR excess emission fr om circumstellar di sks which could 
bias the photometry (jStolte et al.|| 2010 |). Adopting t^ 
single power-law model described above (Equations 
[^, we find the power-law slope of stars brighter than 
E153M = 17 mag (M >^13 M©) to be notably steeper 
than the slope of stars between F153M = 17-20 mag 
(^2.5 Mq< M <^13 Mq). a Kolmogorov-Smirnov test 
finds the probability of these profiles being drawn from 
the same parent distribution to be <0.05%, demonstrat¬ 
ing that mass segregation is present throughout the spa¬ 
tial extent of the cluster. The profiles in these different 
magnitude bins are shown in Figure with the fit sum¬ 
marized in Table and accompanyin^ivariate posterior 
distributions in Appendix [C] 
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Table 4 

EFF87 Profile Fit Results 


Parameter^ 

Prior^ 

Result 

7 

U(0.5, 4.5) 

2.3 ± 0.2 

b (stars pc“^) 

U(0, 8) 

3.3 ± 1.2 

a (pc) 

— 

0.13 ± 0.03 

So (stars pc“ 

- 

2209 ± 929 

^ Description 

of parameters: 

7 = Outer 


power-law slope; b = field contamination, a 
= core radius, Sq = Normalization factor 
^ Uniform distributions presented as U(min, 
max), where min and max are bounds of the 
distribution 




Figure 11. The radial profile of the Arches Cluster. Left: The best-fit power law model as described in the text and Table The red 
line represents the power-law fit to the unbinned data and the black dotted line the residual field contamination value. A binn^ profile is 
included to guide the eye; the black open points are the binned profile before completeness correction and the red solid points the binned 
profile after completeness correction. Uncertainty in the fit (Icr) is captured by the red shaded region, which spans the standard deviation 
of 1000 profiles randomly drawn from the joint posterior distribution. Note that the binned pro f ile is presented only for comparison and 
does not affect the fit. Right: Same as left, but with the inner cluster profile fromiFspinoza et al.|(|2009|> added in bl ue (10 Mq < M < 120 
Mq) and scaled to our profile. The radius range probed the previous proper motion study ot [uiarkson et al.| ( |2012| > is shaded in grey. 


The adopted magnitude separation is an optimization 
between obtaining a large enough sample for good statis¬ 
tics in the bright-star profile and showing the mass seg¬ 
regation, which become less evident with fainter magni¬ 
tude cuts. To demonstrate this, we split the sample into 
three subsets by magnitude such that each magnitude 
bin contains ^130 cluster members before completeness 
corrections in Figure These magnitude bins corre¬ 
spond to F153M < 17.3 mag (M >^12 Mq), F153M = 
17.3 - 18.8 mag (-6 Mo< M <-12 Mq), and F153M 
= 18.8 - 20 mag (—2.5 Mq< M <—6 Mq). Mass segre¬ 
gation remains evident with the brightest (thus highest 
mass) profile being noticeably steeper than the other two 
profiles, while the intermediate and faint-star profiles are 
more similar to one another. 

We caution that the conversion from observed mag¬ 
nitude to mass is highly uncertain due to uncertainties 
in evolutionary models, especially in the pre-main se¬ 
quence. To determine the stellar masses at the magni¬ 
tudes presented above, we adopt a cluster isochrone with 
the nominal properties of the Arches cluster (age = 2.5 


Myr, distance = 8000 pc, = 2.4 mag) constructed 
using a combination of Ge neva models with initial rota¬ 
tion speed of 300 km s“^ ( Meynet fc Maedei]|2003 ) and 
Si ess et al. ( 2000 ^ pre- main sequence models, as discussed 


Lu et al. 




However, a more accurate conversion 
from magnitude to mass (and a more detailed examina¬ 
tion of cluster mass segregation) will be the focus of a 
future paper. 

3 . 2 . 2 . The Search for Tidal Tails 

Given the strong gravitational fields near the GG, the 
Arches cluster is expected to have tidal tails leading and 
trailing its orbit. Such structures have been observed 
for globular clusters and have yielded insight to the ob- 


(e.g., lOdenkirchen et al. 2001 

2003 Grillmair & John- 

son 2006|). Adopting the mode] 

ot the initial conditions ot 


Glarkson et al. (|2012|), and assuming a current position 
100 pc in front of the GG, iHabibi et al.l (|2014|) predict 


that the Arches should have tidal tails extending 20 pc 
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Figure 12. Mass segregation in the Arches cluster. Top: The clearest evidence for mass segregation is found in the steepening of the 
profile for high-mass stars (black; F153M <17 mag, or M >~13 Mq) compared to low-mass stars (red; F153M >17 mag, or M <^^13 Mq). 
On the left^ the solid lines show the power law profile fit to the unbinned data, the dotted lines show the residual field contamination values, 
and the data points show the binned profiles after completeness correction. Icr model uncertainties are shown as the shaded regions. The 
posterior distributions of the power-law slope F are shown to the right. The slopes differ by 2.5cr, and a Kolmogorov-Smirnov test rejects 
the hypothesis that the profiles are drawn from the same parent distribution. Bottom: Similar to above, but splitting the sample into three 
subsamples by magnitude such that each magnitude bin contains ~130 stars before completeness correction. Mass segregation remains 
evident in the steepening of the brightest profile compared to the other two profiles. The best parameter values for all fits are presented in 
Table O 


(^500”) along the Galactic plane. To search for these 
structures we compare the radial profiles parallel and 
perpendicular to the cluster’s bulk velocity, whi ch is con¬ 


sistent with the direction of the Galactic plane (Glarkson 


et al. 


2012 ). 


T'idal tails would cause an asymmetry in these profiles, 
either as a steepening or truncation of the perpendicu¬ 
lar profile relative to the parallel profile as the cluster 
is stretched and sheared by the Galactic tidal field (i.e. 
Figure 5 of jOdenkirchen et al.| 2003). Using the single 
power-law rnodel, we do not hnd evidence of a significant 


difference between the parallel and perpendicular profile 
slopes (Table Figure [T^. Bivariate posterior distribu¬ 
tions for these profiles are are presented in Appendix 0 
While the binned profiles might appear to be discrepant 
at - 1.5 pc, a Kolmogorov-Smirnov test of these pro¬ 
files in the region of highest completeness (0.5 pc < R 
< 3.0 pc) concludes that the profiles are drawn from the 
same distribution with a probability of ^16%. Thus, we 
cannot conclude that profiles are statistically different. 
This conclusion does not change when different cluster 
membership probability cuts are adopted (P > 0.7, for 
example). Further observations of the Arches cluster at 


large radii are needed to detect the presence of tidal tails. 


3.3. Observed Tidal Radius 

Utilizing the large field of view, we can directly con¬ 
strain the spatial extent of the Arches cluster for the 
first time. This is a significant improvement over previ¬ 
ous studies which were forced to estimate the tidal radius 
based on obs ervations o f the in ner region of the cluster. 
For example, Kim et al. (2000) compared the radial pro¬ 
file of massive stars m the Arches (M > 20 Mq) out to 
0.8 pc (Figer et al. 1999a) to the radial profiles of simi¬ 
larly massive stars m M-body simulations of the cluster 
and found an expec ted tidal radius between 1-1.2 pc. A 


second estimate by Portegies Zwart et al. (2002) placed 
the tidal radius at 1.6 - 2.5 pc, based on a highly model- 
depende nt analysis of the m ass segregation observed in 
the same Figer et al. (1999a) profile. However, our study 
shows that the Arches profile extends well beyond these 
predictions, with no evidence of King-like tidal radius 
out to ^3 pc (Figure p^. 

To place a quantitative lower limit on the observed 
tidal radius, we use the Bay esian framew ork described 
above to fit our profile with a King (1962) model: 
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Figure 13. The search for tidal tails in the Arches cluster. Left: No significant asymmetries suggesting the presence of tidal tails are found 
in a comparison of the profile parallel (red) and perpendicular (bla ck) to the cluster’s orbit. Power-law fits, residual field contamination, 
and Icr uncertainties are shown in the same manner as Figure [TT] while the best-fit parameter values are presented in Table Right: 
The posterior distributions for the power-law slope T of the fitted profiles. The black and red dotted lines show the best-fit slopes for the 
perpendicular and parallel profiles, respectively, and which only differ by ~ 1.4(7. A Kolmogorov-Smirnov test cannot reject the hypothesis 
that these profiles were drawn from the same parent distribution. 


S(r) = k ^ 


([l + (r/rc)2]i/2 


1 

[l + (n/re)2]l/2 



where /c is a normalization constant, 6 is a constant back¬ 
ground term, and Vc and Vt are the core and tidal radii of 
the cluster, respectively. We adopt the core rad ius of 0.14 
± 0.05 pc measured by Espinoza et al. (2009) as a prior 
for Tc, though our profile provides no additional infor¬ 
mation at R < 0.25 pc. An uninformed prior is used for 
Vt , and k is calculated such that the integral of the fitted 
profile yields the total number of cluster members ob¬ 
served after membership probability, completeness, and 
area corrections. 

The result of the fit is summarized in Table |5] and 
marginalized posterior distributions presented in Ap¬ 
pendix 0 We obtain a 3cr lower limit of 2.8 pc on a 
King-like tidal radius for the cluster. Of course, it is 
quite possible that the Arches profile is truncated at all, 
and may behave as a power law throughout the full clus¬ 
ter extent. Regardless, it is clear that the cluster extends 
beyond its predicted tidal radius of 2.5 pc. 


3.4. Tidal Breaks in the Radial Profile 

We examine our profile for the presence of breaks which 
would be indicative of significant tidal interactions. The 
profile is remarkably consistent with a single power law 
from 0.25 pc < R < 1 pc, though beyond 1 pc the pro¬ 
file appears to exhibit slightly higher stellar densities 
than expected. This feature is adequately modeled as a 
constant field contamination term in Equation How¬ 
ever, we must confirm that this feature can indeed be 
attributed to residual field contamination rather than a 
true over-density of cluster stars. 

Eield contamination may arise from the uncertainties 



Figure 14. The background-subtracted profile of the Arches clus¬ 
ter compared to a King profile with a tidal radius of it = 2.5 pc. 
The Arches profile is composed of the best-fit power law slope (red 
line) and binned profile (black points ) from this study and t he pro¬ 
file for the inner 0.4 pc measured by|Espinoza et al.|(|2009| scaled 
to our profile (blue line). The King prohle model (black dotted 
line) is clearly discrepant with this profile at large radii. We place 
a 3(7 lower limit of 2.8 pc on the location of a King-like tidal radius 
in the Arches cluster. 


in the cluster membership probabilities, which in turn 
are a result of uncertainties in the fits of the cluster and 
field kinematic distributions. To estimate the impact of 
these uncertainties, we perform a Monte-Carlo experi- 
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Table 5 

King Profile Fit Results 


Parameter^ 

Priori 

Result 

n (pc) 

U(l, 15) 

2.8^ 

Tc (pc) 

G(0.14, 0.05)'^ 

0.13 ± 0.03 

b (stars pc“^) 

U(0,15) 

3.46 ± 0.94 

k (stars pc“^) 

— 

1729 ± 643 


^ Description of parameters: rt = tidal radius; 
Tc = core radius; b = field contamination; k 
Normalization factor 

^ Uniform distributions: U(min, max), where 
min and max are bounds of the distribution; 
Gaussian distributions: G(mean, standard de¬ 
viation) 

3f7 lowei ^ limit 
^ Source: 


Espinoza et al. (2009) 


merit where we randomly draw 1000 kinematic mo dels 
from the joint posterior distribution obtained in § |2.2| and 
calculate the stellar cluster membership probabilities for 
each. We then determine the number of cluster members 
for each model from the sum of membership probabilities 
for all stars with Pmember >0.3, the sam e cal culation we 
do for the best-fit kinematic model in § |3.1[ The stan¬ 
dard deviation in the number of cluster members across 
the kinematic models is 34.32, centered around a median 
very nearly equal to the number of members identified by 
the best-fit kinematic model. We therefore adopt 34.32 
as the uncertainty in the number of cluster members due 
to imperfect cluster membership probabilities. 

The cluster membership uncertainty provides an es¬ 
timate on the number of field contaminants potentially 
among our sample of cluster members. Spread evenly 
across the field, this would result in a field surface den¬ 
sity of 1.49 stars pc“^ in the profile. This is consistent 
with the background of the single po wer law fit (2.52 T 
1.32) to Icr. In addition. Figure 1 of Peharrubia et al. 
(2009) shows that a tidal event may cause an asymme¬ 


try m the parallel and perpendicular profiles, with the 
stellar over-density dominating the parallel profile rela¬ 
tive to the pe rpend icular one. That no such asymmetry 
is observed ( p.2.2 ) is further evidence that there is no 
tidal break in the Arches profile between 0.25 - 3.0 pc. 

4. DISGUSSION 

In order to interpret the Arches cluster profile in terms 
of its tidal history, we require theoretical/numerical stud¬ 
ies of YMCs on moderately eccentric orbits in the inner 
Milky Way potential that examine the evolution of the 
outer radial profile. If it exists, the observable King-like 
tidal radius of the cluster is significantly larger than pre¬ 
dicted by previous studies, though these assume a spher¬ 
ically symmetric potential for the Galaxy or assume the 


cluster is on a circ ular orbit (Kim et al 


Zwart et al.||2002 ). Given the high gas densities toward 


2000 Portegies 


the GC, the ettect of interactions with Giant Molecular 
Glouds on the radial profile may also need to be consid¬ 
ered. Unfortunately there are no such studies currently 
available in the literature; while there are many N-body 
simulations of clusters in varying tidal fields, very few 
discuss the correspondi ng e volution of the outer cluster 
profile. As a result, in ^4.1|we consider the implic ations 


of simulation results by Peharrubia et al. (2009 here- 


galaxies to tidal perturbations at perigalacticon. We 
discuss other simulations which more closely reflect the 
Arches cluster and its environment but do not examine 
the e volution of the radial profile in sufficient detail in 

El 

4.1. Time Since Last Pericenter Passage 

Modeling the effects of tidal stripping on dwarf 
spheroidal galaxies, P09 find a relation between the lo¬ 
cation of a tidal break in a radial profile and the time 
elapsed since perigalacticon. Here we apply this rela¬ 
tion to the Arches cluster to place limits on the time 
elapsed since its closest approach to the GG. However, 
there is a major caveat in this analysis: the orbits ex¬ 
amined by P09 are significantly different than what is 
expected for the Arches. These authors simulate dwarf 
spheroidal galaxies on highly eccentric orbits (e = 0.96 
- 0.99) with a closest approach of 900 pc from the GG, 
where the Arches is likely on an orbit with e = 0.25 - 
0.38 and a closest approach between 50 - 200 pc from the 
GG ( Stolte et al.||2QQ8 Krnijssen et al.||20I5 ^ As sumin g 
the Milky Way potential model described in Bovy (2015), 
the radial force felt by the Arches is ^6 times larger than 
that of the innermost POO orbit. A similar formation and 
evolution of a tidal break in response to perigalacticon is 
found in N-body simulations of dwarf galaxies o n more 
moderately -eccentric orbits (0.23 < e < 0.9) by Lokas 
et al. (20I3|), though these models also probe weaWr tidal 
helds than is experienced by the Arches (smallest peri¬ 
center: 12.5 kpc). N-body studies of star clusters on 
eccentric orbits show the formation of a tidal break at 


1999 


Lee et al. 


large radii, as well (Khpper et al. 2010 Johnston et al 


2005r Given the supporting evidence 


from multiple studies and the lack of alternative studies 
in stron ger ti dal fields that examine the outer radial pro¬ 
file (see p.2[ ) , we move forward with the results from the 
POO simulations. 

After perigalacticon, the radial profiles of POO develop 
a tidal break that initially forms at small radii and moves 
outward over time. Normalizing by core radius Rc and 
core crossing time ter = Rc / they obtain the follow¬ 
ing relation: 


Rb/Rc = 0.55(i - tp)/tc 


( 10 ) 


after POO), which model the response of dwarf spheroidal where Rf, is the radius of the tidal break and (t — tp) 
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is the time elapsed since perigalacticon. Given that no 
break is observed in the Arches profile between 0.25 pc - 
3.0 pc and adopting a velocity dispersion of 5.4 km s~^ 
(Clarkson et al.|2012), this relation indicates that the 
Arches has not had a significant tidal perturbation be¬ 
tween ^0.08 Myr and Myr ago. However, POO note 
that while their relation describes all models well at large 
break radii, it tends to underestimate the time since peri¬ 
galacticon passage for break radii close to the core radius 
( Rb < 4 Rc). Restricting the lower boundary to the 
innermost radius that is well described by this relation 
(0.6 pc for the Arches), we find that it is likely that the 
Arches has not experienced perigalacticon between ^0.2 
Myr and Myr ago. 

This result, which suggests that the cluster may ei¬ 
ther be nearly at or long past closest approach, places 
a limit on the set of potential orb its calculated for the 


Arches cluster (Stolte et al. 
bits to those which place t 


[2008). Restricting the or- 
ie Arches within the Cen¬ 


tral Molecular Zone (highly likely given its inte ractions 
with surrounding gas clouds, cf. Wang et al. 2010), 


this constraint rejects the viability of retrograde orbits 
since these place the cluster’s closest approach within the 
timeframe in which we would expect to observe a tidal 
break. This provides further evidence that the cluster is 
on a prograde orbit and thus is located in front of the 
sky plane which passes through Sgr A*. However, we 
cannot significantly constrain the prograde orbits, which 
place the cluster very near closest approach. The al¬ 
lowed orbits are consistent with both cluster formation 
mechanisms discussed in ^ and so no additional insight 
regarding the birth of the Wches cluster can be obtained. 

Additional astrometric observations of the Arches are 
needed to provide higher precision proper motions from 
which the velocity dispersion profile of the cluster can be 
measured. Combined with the radial profile, the velocity 
dispersion profile is a sensitive tracer of cluster’s tidal in¬ 
teraction history and current dynamical state, and may 
lead to the measurement of current Jacobi radius of the 


cluster ( jKupper et aL][^lQ| ). This can yield the present 
distance between the Arches and Sgr A*, the last bit of 
information required for a full orbital solution. The dif¬ 
ferent formation mechanisms for the Arches can then be 
distinguished, as they differ in predictions of the current 
distance between the cluster and SMBH (^50 pc for the 
tidal compression scenario versus ^100 - 200 pc for the 
cloud collision scenario). 

It is important to note that the orbit calculations of 
the Arches cluster rely on the accurate measurement of 
the cluster’s bulk motion relative to Sgr A*. This is not 
a trivial task due to an observational bias towards stars 
on the near side of the GC. Limited by their field of view, 
Stolte et al. (2008) model the kinematics of 67 identified 
held stars using a single circular Gaussian distribution, 
taking the bulk motion of the Arches to be the difference 
between the kinematic centers of the Arches and field 
distributions (5.6 ± 0.5 mas yr“^). This is the value 
used in the currently published orbits. Using 210 field 
stars and a more sophist icated elliptical Gauss ian model 
for the field population, Glarkson et al. (2012) obtain a 
slightly lower field motion of 4.39 ± 0.38 mas yr“^ rel¬ 
ative to the cluster. However, with a much larger field 
of view and more available field stars (^5322 stars with 


Fmember < 0-3), our study shows that the field has a 
complex kinematic structure that must be modeled with 
multiple Gaussians. Great care must be taken to prop¬ 
erly interpret these structures in the context of a Galactic 
model and measure the bulk motion of the Arches with 
respect to Sgr A*, which is left to a future paper. 

4.2. Applicability of Other Theoretical Studies 

Despite the significant caveats in applying the POO re¬ 
sults to the Arches cluster, it is the most applicable the¬ 
oretical study currently available that examines the de¬ 
tailed evolution of the outer radial profile of a stellar 
system as it passes through perigalacticon. N-body sim¬ 
ulations of dwarf galaxies on or bits with moderate eccen¬ 
tricities by Lokas et al. (2013) exhibit tidal breaks that 
behave similarly as those in the POO models, though the 
relation between the location of the break radius and the 
time since perigalacticon passage is not assesse d quan- 
titively. Al ternatively, N-body simulations by Kfipperj 
et al. (2010) examine objects more similar to the Arches, 
studying the behavior of 10 ^ Mq star clusters on elliptical 
orbits with 0.25 < e < 0.70 (smallest pericenter: 600 pc). 
The resulting cluster profiles show that extratidal stars 
form power-law extensions at large cluster radii, though 
the detailed evolution of the profile after perigalacticon 
is not explored in detail. Additional simulation s of clus¬ 


ters on ellipt ical orbits by Johnston et al. (1999) and Lee 
et al. (2006) also show extratidal stars forming a break 


m the cluster profile, though neither study examines how 
this break evolves as a function of orbital phase. 

There is a large body of additional literature studying 
the dynamical evolution of star clusters in tidal fields, 
though these do not show the evolution of the outer ra¬ 
dial profil e. N-body simulations of clust e rs on eccentric 


orbits by B aumgardt fc Makin o (2003), Lamers et al. 
(2010), and" Webb et al. ( [2014 ) primarily focus on the 


evolution of the m ass-loss rate and mass function, while 
Webb et al. (2013) examines the half-mass radius, tidal 
radius, and cluster size rather than the morphology of the 
radial profile at large radii. Many other studies exam¬ 
ine clusters on circular orbits, th ough these are liinited to 
old (>10 G yr) globular clusters (Trent! et al.|2010 [Gieles 


et al.||2011|), or do not present detailed radial profil es of 
their models ([Ernst et al.||2009t [Madrid et al.[[2012 ). 

It is worthwhile to note that none of the studies dis 


cussed above examine the evolution clusters within the 
central regions of the Milky Way. There have been sev¬ 
eral studies of Arches-like young compact clusters within 
the inner 200 pc of the Galaxy, but these similarly do not 
examine the ev olution of th e outer radia l profile. N-body 


simulations by Kim et al. (2000) and Portegies Zwart 
et al.[ ([2002[) make predictions regarding expected clus¬ 
ter lifetimes and the evol ution of the radial prof i le out 
to the half-mass radius. Portegies Zwart et al. (2004) 
focus on the evolution of the mass function of Arches- 
like clusters, concluding that the ma ss function of the 
inner region of the Arches reported by Eiger et al. (2002) 
could be explained by dynamical mass segregation, bev- 
eral other studies model the effects of dynamical friction 
on compact clusters, predicting in-spiral towards the GG 
and their subsequent evolution, though th e cluster pro¬ 
file during this process is not presented ([Kim fc Mor- 
ns||2QQ3l [Portegies Zwart et al.[[200^ IGiirkan ^ Kasio 
2005p . Additional studies of the behavior of the radial 
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profiles of Arches-like clusters near the GC are required 
to draw more conclusive interpretations from the obser¬ 
vations presented here. 

5. CONCLUSIONS 

We have conducted a multi-epoch photometric and as¬ 
trometric study of the Arches cluster using the Hubble 
Space Telescope WFC3IR camera at 1.27, 1.39, and 1.53 
/im. Using a sophisticated astrometric pipeline we ex¬ 
tract individual stellar proper motions to an accuracy of 
at least 0.65 mas yr“^ down to F153M ^ 20 mag (^2.5©), 
reaching a precision of ^0.1 mas yr“^ for the bright¬ 
est stars. Taking advantage of the distinct kinematic 
properties of the cluster, we use a 4-Gaussian mixture 
model to simultaneously fit the cluster and field proper 
motion distributions and calculate cluster membership 
probabilities. This is a substantial improvement over 
photometrically-determined cluster membership due to 
the large degree of differential reddening across the field. 
The field of view in this study is 144 times larger than 
previous astrometric studies of the Arches cluster, al¬ 
lowing for the identification of high-probability cluster 
members out to a cluster radius of 75” (^3 pc at 8 kpc). 

Gombining the cluster membership probabilities, an 
extinction map derived from red clump (RG) stars, and 
an extensive completeness analysis, we construct the stel¬ 
lar radial density profile for the Arches cluster between 
6.25” < R < 75” (0.25 pc < R < 3.0 pc) down to a dif¬ 
ferentially de-reddened magnitude of F153M = 20 mag. 
This profile is well fit by a single power-law of slope F = 
2.06 ± 0.17 with a constant field contamination density 
of 2.52 ± 1.32 stars pc“^. Surprisingly, no evidence of 
a tidal radius is observed. Adopting a King profile as a 
model, we obtain a 3cr lower limit of 2.8 pc for the ob¬ 
served tidal radius of the Arches cluster. This shows that 
the cluster extends beyond its largest predicted theoret¬ 
ical tidal radius of 2.5 pc. 

Additionally, we examine the Arches cluster profile for 
evidence of mass segregation and tidal tails. We find the 
cluster to exhibit mass segregation at all observed radii, 
with the radial profile of bright stars (F153M < 17 mag, 
or M> ^13 Mq) being notably steeper than the profile of 
fainter (17 < F153M < 20 mag) stars. A KS test reveals 
the differences between these profiles to be significant. 
We leave a careful conversion from brightness to mass for 
a future paper. Similarly, we search for evidence of tidal 
tails by comparing the profile parallel to the direction of 
orbit to the profile perpendicular to it. No statistically 
significant asymmetries are observed in these profiles, as 
would be expected from tidal tail structures. Further 
observations, perhaps at larger cluster radii, are needed 
to continue to search for tidal tails. 

No evidence of a tidal break is observed in the ra¬ 
dial profile, as might be expected if the Arches has ex¬ 
perienced a tidal perturbation in its recent past. As¬ 
suming that the results of dynamical simulations of 
dwarf spheroidal galaxie s on highly eccentric orbits by 
Peharrubia et al. (2009) can be applied to the Arches, 
this suggests that the Arches not likely experienced its 
closest approach to the GG within 0.2 - 1 Myr ago. This 
constraint would reject all possible retrograde orbits of 
the cluster, providing further evidence that the Arches 
is on a prograde orbit and located in front of the sky 
plane which intersects Sgr A*. However, additional sim¬ 


ulations studying the profile of Arches-like clusters on 
mildly-eccentric orbits in the inner Milky Way potential 
are required to interpret the observed profile with higher 
confidence. 

Further astrometric observations of the Arches to ob¬ 
tain its velocity dispersion profile are needed to better 
constrain its orbit and distinguish between different pos¬ 
sible cluster formation scenarios. It is important to note 
that an accurate determination of the cluster’s orbit re¬ 
quires measuring its bulk motion with respect to Sgr 
A*, which is a difficult task given the complex kinematic 
structure of the field population revealed in this study. A 
revised measurement of the Arches proper motion, along 
with a new calculation of possible orbits, is left to a fu¬ 
ture paper. 
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APPENDIX 

A. WFC3IR MEASUREMENTS AND PROPER MOTIONS 

In this appendix we describe the methods and software used to extract high precision astrometry, photometry, and 
proper motions from the WFC3IR observations. This is the first application of this methodology on WFC3IR observa¬ 
tions. Stars are first detected and measured in the fit images using the FORTRAN p rogram imq2xym-wfc3ir stdpsf 
developed by Jay Anderson. Similar to the code img2xym-WFC. 09x10 documented in Anderson & King (2006), this 
performs PSF-fitting measurements using a library of spatially-variable PSF models it derives tor the Wl'(J3iR camera. 
This library contains a 3x3 grid of PSFs that spans the camera’s field, where the PSF at any point can be derived 
from a spatial interpolation of these models. Since this routine operates on one image at a time in a single pass, it 
is not designed to deal with overlapping stars. It reduces each star as if it is the only contribution to the 5x5 pixels 
centered on its brightest pixel. 

Each image is run through this program twice. The first iteration extracts the bright high S/N stars in the field 
(^400 stars per image), simultaneously measuring the residuals between the observed PSF and the library PSF. New 
image-specific PSF models are created in order to minimize these residuals and make them uniform across the field. 
The second iteration then uses the modified PSF library to accurately measure both bright and faint sources, producing 
a star list with fluxes and positions for ^13,000 sources in each image extending down to FI27M = 22.45 mag, FI39M 
= 22.09 mag, and FI53M = 21.71 mag. A small number of stars (^200) with FI53M <^I5 are saturated and are 
thus measured using the outer part of the PSF. 

Next, we cross-identify stars from each exposure with a master list for the filter/epoch set, initially taken to be 
the first image in the set. Common stars found in at least 75% of the images are used to transform positions from 
the distortion-corrected frame of each exposure into the master list reference frame using general 6-parameter linear 
transformations. This gives ^N observations for each star, where N is the number of images in the filter/epoch 
set, which allows us to find a robust average position and flux. These averaged measurements produce an overall 
star catalog for the filter/epoch. We then adopt these catalogs as the new reference-frame positions and repeat 
the procedure, improving the transformations. The improvement in the second iteration is considerable, decreasing 
astrometric residuals by nearly a factor of 2. Finally, the new star catalogs for each filter/epoch are then transformed 
to an arbitrary astrometric reference frame where the net motion of the cluster plus field stars is 0 mas yr“^. This 
produces what we will call the “one-pass” catalogs, because they are limited to the stars which can be detected in a 
single image. 

In principle, we could continue our analysis of the Arches Cluster using the one-pass catalogs. However, by stacking 
the images in each filter/epoch we significantly increase the detection depth, important because signal from otherwise 
undetected faint stars can be mistakenly associated with brighter stars and introduce biases in the measurements. 
With the image transformations, PSF models, and one-pass catalogs as input, we use the program KS2 to stack the 
images, make stellar detections to significantly fainter magnitudes, and then redo the astrometric and photometric 
measurements in each individual image at the positions found in the stacked image. The measurements for the 
individual images are then averaged together to produce the final catalogs for each epoch/filter. We emphasize that 
we don’t use the measurements of the stacked image, but only those of the individual images. This iterative multiple¬ 
finding strategy increases the number of stars detected to ^50,000 per filter (nearly five times as many as th e one-pass 
analysis), reaching FI27M = 23.63 mag, FI39M = 23.29 mag, and FI53M = 23.31 mag as noted in 


2.1 


To ensure 

the accuracy of our transformation we again combine the individual star lists into star catalogs for each falter/e poch 


and re-align to the common astrometric reference frame using a first-order bi-variate polynomial (see Ghez et al. 2005 
and references therein). Higher order fits were found to introduce artificial structure in the astrometry. Astrometric 
and photometric errors as a function of magnitude are presented in Figure 

The KS2 code returns a final star catalog with the root-mean-square errors {(Trms) for the astrometric and photo¬ 
metric measurements for each star in the filter/epoch. Theoretically the errors should be quantified by the error on 
the mean {(Trms / V^obs , where Nobs is the number of images in the filter/epoch). To test this, we compare the KS2 
errors to quantities measured across the FI53M epochs that are directly caused by these errors: the standard deviation 
of the magnitude of each star (assuming no stellar variability) and the residuals between the measured position and 
position predicted by the star’s proper motion. This comparison reveals the error on the mean to better capture the 
astrometric errors and the RMS error to better capture the photometric errors, and so we adopt these for the individual 
star measurements throughout (Figure[^. This choice does not significantly affect the proper motion errors described 
below, which are dominated by the position residuals rather than the astrometric errors themselves. 

Proper motions are derived independently for the X and Y directions (in the image coordinate system) using a linear 
fit to the change in position over the FI53M epochs: 


X = xoFVx{t-to) (AI) 

y = yoFvy{t-to) 


where to is the astrometric error-weighted average time of observations, (v^^Vy) are the X and Y proper motions, 
and (to, ^o) and (t, y) are the star positions at to and t, respectively. To test the validity of these errors, we 
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Figure 15. Astrometric error on the mean and photometric RMS error vs. observed magnitude for the F127M (blue), F139M (green), and 
F153M (red) filters. The solid lines show the median errors and the shaded regions cover one standard deviation. Stars with astrometric 
errors above 1.5 mas (and thus proper motion errors above 0.65 mas yr“^) or photometric errors above 0.06 mag in the F153M filter are 
not included in the analysis. These cuts are shown by the black dotted lines. 


apply the derived proper motions to their respective stars and examine the residuals between the predicted and the 
observed positions. The distributions of the X and Y position residuals are approximately Gaussian, though more 
power is present in the wings than expected (Figure fl^. This is a consequence of stellar crowding distorting our 
measurements, as these wings are dominated by faint stars (F153M < 20 mag) which are more prone to this effect. 
Given the high stellar density of near the center of the cluster this is not unexpected. The distribution of values 
for the proper motion fits follows the expected distribution for 1 degree of freedom, appropriate as we constrain 2 
parameters with 3 measurements. 

Photometry is calibrated to the standard Vega magnitude system by deriving the filter-dependent offset betw een KS2 
magn itudes and 0.4” aperture photometry magnitudes. Aperture photometry is performed using DAOPHOT (Stetson 
|1987) on the drz image for each epoch/filter, a composite image of all exposures in the same filter/epoch produced by 
the n ST p ipeline. This offset is then combined with the appropriate 0.4” zeropoint derived for the WFC3IR camera 
(see for reference) to determine the overall z erop oint for the KS2 observations. For stars with F153M < 20 


(consistent with our proper motion precision cut, § 2.1), the median F153M astrometric and photometric errors are 


0.34 mas and 0.018 mags, respectively, with evidence of higher errors in regions of increased stellar density (Figure 
|T7| ). Within this sample the median F127M and F139 photometric errors are 0.033 mag and 0.025 mag, respectively. 

B. ARTIFICIAL STAR AND OBSERVED STAR ERRORS 


The completeness analysis described in § |2.4| assumes that the measured artificial star errors match the observed star 
errors. A direct comparison of the errors reveals that the observed astrometric and photometric errors have an error 
floor that is not reproduced by the artificial star tests (Fig. |T^ . A possible explanation for this feature is residual 
PSF variations which are not captured by our spatially-varying PSF model, as the artificial stars are planted using 
this model and thus wouldn’t reflect this error. A constant error correction term of 0.14 mas and 0.008 mag is added 
to the artificial star position and magnitude uncertainties, after which the artificial star errors are found to closely 
follow those of the observed stars. This error correction does not have a noticeable effect on the completeness analysis 
as only stars far below the error cuts are significantly affected. 

C. POSTERIORS OF PROFILE FITS 

I n th is appendix we present the bivariate posterior distributions for the different power law profile fits des cribed 
in §3.2| and the marginalized and bivariate posterior distributions for the King profile fit described in §3.3[ The 
bivariate posterior distributions show the correlations between the power law slope (F), background level (b) and 
profile amplitude (Aq) for the profile fits to the full cluster (Figure!^, high-mass/low-mass members (Figure [^, 
and the parallel/perpendicular members (Figure 22). The marginalized posteriors of these parameters (not shown) 
are well described by Gaussians, which provide theoest fit parameter values and errors reported in Table 
For the King model fit, the marginalized posterior distribution for the tidal radius (r^), core radius (rc), background 
(5) , and normalization factor {k) are shown in Figure 23, while the bivariate po steri or distributions are shown in Figure 
24 The significance of the 3cr lower limit on the tidal radius is discussed in 13.3 The best fit parameter values and 
errors are reported in Table 
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Figure 16. Top: The fractional difference between the RMS error {left) and error on the mean {right) and the “true” errors for the 
photometric measurements. The “true” photometric error for each star is taken to be the standard deviation of the observed magnitude 
across the F153M epochs. Bottom: Similar to above, the fractional difference between the RMS error {left) and error on the mean {right) 
and the “true” errors for the astrometric measurements. The “true” errors are taken to be the RMS residuals between the fitted proper 
motion and the observed position of each star in the F153M epochs. For the photometry, the mean is found to underestimate the true 
error, and so the photometric RMS error is adopted for individual measurements. For the astrometry, the RMS error is found to strongly 
overestimate the true error, and so the astrometric error on the mean is adopted for individual measurements. 
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Figure 17. Average astrometric error {left) and F153M photometric error {right) as a function of position on the camera for all stars 
with F153M < 20 mag. Average error values are calculated in 7” bins and are plotted relative to the cluster center. Axes are oriented in 
the same manner as Figure Higher errors are observed in the low reddening regions and the cluster center due to stellar crowding. 
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Figure 18. Positions residuals between the observed positions and those predicted from the fitted proper motions. Top: A histogram 
of the X (left) and Y (right) residuals (blue line). Ideally these distributions would be Gaussian (red line), but stellar crowding results 
in more power in the wings of the distribution. Bottom: The distribution of values for the proper motion fits in both X (left) and Y 
(right)^ where we adopt the error on the mean as the proper motion error (blue line). These values match the expected x^ distribution 
with 1 degree of freedom (red line), validating our reported errors. 
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Figure 19. Comparison of the astrometric (left) and photometric (right) errors for the observed and artificial stars in blue and red, 
respectively. The top plots show the initial artificial star errors and the bottom plots show the artificial star errors after a constant error 
term has been added. The solid line and filled area represents the median and standard deviation of the errors. The additional error terms 
refiect the error fioor seen in the observed star measurements, and are measured to be of 0.14 mas and 0.008 mag for the astrometry and 
photometry, respectively. After this adjustment the error distributions match well. The green dotted line shows the error cuts used in this 
study. 
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Figure 20. The bivariate posterior distributions for the power law model fit to the full cluster sample. 
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Figure 21. The bivariate posterior distributions for the power law model fits to the high-mass (top, F153M <17 mag) and low-mass 
{bottom, F153M > 17 mag) cluster members. That the high-mass stars have a steeper power law slope F is an indication of mass segregation 
in the Arches cluster. 
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Figure 22. The bivariate posterior distributions for the power law model fits to the cluster members parallel (top) and perpendicular 
(bottom) to the bulk cluster orbit. The existence of tidal tails would cause asymmetries in these profiles such as a difference in the power 
law slope F of these profiles. No significant evidence for tidal tails are found. 
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Figure 23. The ID marginalized posterior distributions for the King model fit. The output of the Multinest sampling is in black, the 
corresponding Gaussian fit in red, and the input prior in green (if applicable). The 3a limit to the tidal radius rt (2.8 pc) is indicated by 
the red line, where 99.7% of the best-fit models fall above this value. 



k (stars pc k (stars pc b (stars pc 


Figure 24. The bivariate posterior distributions for the King model fit. 




























